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Abstract 

An exactly-solvable model for the decay of a metastable state coupled to a semi-infinite tight- 
binding lattice, showing large deviations from exponential decay in the strong coupling regime, is 
presented. An optical realization of the lattice model, based on discrete diffraction in a semi-infinite 
array of tunneling-coupled optical waveguides, is proposed to test non-exponential decay and for 
the observation of an optical analog of the quantum Zeno effect. 
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The understanding and control of the decay process of an unstable quantum state has long 
been a subject of debate in different areas of physics. Though an exponential law is known to 
be a good phenomenological fit to many decay phenomena, quantum mechanics ensures that 
;he survival probability P(t) is definitely not exponential at short and long times (see, e.g., 
M, l2l, [3j). In particular, at short times P(t) always shows a parabolic decay, i.e. dP/dt — > 
as t — > 0. These universal features have been extensively investigated in some specific models 



% [5| , or in the 
which describes 



describing the tunneling escape of a particle through a potential barrier 
framework of the exactly- solvable Friedrichs-Lee Hamiltonian jfj, R S. R [lO ] . 
the decay of a discrete state coupled to a continuum. The short-time features of the decay 
process have attracted much attention because they can lead, under certain conditions, to 
either the deceleration (Zeno effect) or the acceleration (anti-Zeno effect) of the decay by 
frequent observations of the system (see, e.g., (si, [9, 11] and references therein). Evidences of 



non-exponential decay features at short times and the observation of the related Zeno and 
anti-Zeno effects have been reported in recent experiments on quantum tunneling of trapped 
sodium atoms in accelerating optical lattices [3]. Similar effects have been proposed to occur 



for quantum tunneling in analogous macroscopic systems, such as Josephson junctions [13 ]. 
In this Letter a novel and exactly-solvable model of non-exponential decay of an unstable 
state tunneling-coupled to a tight-binding lattice is presented. A simple and experimentally 
accessible realization of the model, based on discrete diffraction of photons in an array of 
optical waveguides fl^ ]. is proposed along with an optical analog of the quantum Zeno effect. 
To set our model in a general context, we consider a semi-infinite lattice described by the 
tight-binding Hamiltonian [Fig. 1(a)]: 

00 

H TB = -ft £ A n (\n)(n + l\ + \n + l)(n|) , (1) 

n=l 

where \n) (n > 1) is the state localized at the n-th site of the lattice and A n is the hopping 
amplitude between adjacent sites \n) and \n + 1). We assume that for n > 2 the lattice 
is periodic so that, after a rescaling of time t, we may assume A n = 1 for n > 2. The 
boundary site |1) is then coupled to the periodic lattice by a hopping amplitude A x = A, 
which is assumed to be smaller than 1. The tight-binding Hamiltonian (1) has been often 
used as a simple model to describe coherent transport properties and tunneling phenomena 



in different physical systems, including semiconductor superlattices [15(, arrays of coupled 
quantum dots [if]], Bose-Einstein condensates in optical lattices [l^], and arrays of optical 
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FIG. 1: (a) The semi-infinite tight-binding lattice model, (b) Optical realization of the tight- 
binding model based on an array of coupled optical waveguides, (c) Refractive index profile n(x)—n s 
of the waveguide array used in the numerical simulations (parameter values are: n s = 2.138, 
A = 1.55 am, An = 2.4 x 10~ 3 , and a = 12 am). 



waveguides [14] . Il8|. In particular, model (1) can be derived from the continuous Schrodinger 
equation 



n 2 d 2 ^ 



+ V(x)tp, 



(2) 



2m dx 2 

with a potential V(x) = S^Lo Ko(^~ x n ) describing a semi-infinite chain of identical symmet- 
ric quantum wells V w (x) \V W (— x) = V w (x) and V w (x) — > for x — > oo], placed at distances 
%n+i — x„ = a for n > 2 and x 2 — X\ = a > a [see Fig. 1(a)]. If the individual potential well 
V w (x) supports a single bounded mode <p(x) of energy E and if tunneling-induced coupling 
of adjacent wells is weak, Eq.(2) can be reduced to the discrete model (1) by means of a 

nn . n 

tight-binding [17|, [19[ or a variational [20( analysis. After expanding the state of the 
system as \ip) = J2n c nif) exp(— iEt/h)\n), where |n) = <p(x — x n ) is the localized state at 
the n-th well in the chain, in the nearest-neighbor approximation from Eq. (2) one can derive 
the following equations of motion for c n : 



ic-i = — Ac 2 , ic 2 = — c 3 — Aci . 
ic n = -(c n+ i + c n _i) for n > 3, 



(3) 



where A ~ [/ dx <p(x — ao)V w (x)(p(x)]/[f dx <p(x — a)V w (x)<p(x)] is the normalized hopping 
amplitude between states |1) and |2). For A = 0, i.e. for ao/a — > oo, the site |1) is decoupled 
from the other lattice sites and if the system is initially prepared in state |1) it does not decay; 



as A is increased, tunneling escape is allowed and state |1) becomes metastable. The limits 
A — ► and A — > 1 correspond to the weak and strong coupling regimes, respectively. The 
occupation probability of site |1) at time t is given by P(t) = |ci(t)| 2 . Following Gamow's 
approach to quantum tunneling decay [4j, the 'natural' decay rate 70 of state |1), which 
would correspond to an exponential decay law P(t) = exp(—j t), can be readily calculated 
by looking for complex energy eigenfunctions of Htb with outgoing boundary conditions 
(Gamow's states), yielding: 

7o = 2A 2 (l- A 2 )" 1 / 2 . (4) 

However, the exponential decay law turns out to be incorrect, especially in the strong cou- 
pling regime A — > 1 where it fails to reproduce the exact decay law at any time scale. 
According to Ref.Q], one can introduce an effective decay rate 7 e //(t) by the relation 
J e ff{t) = — (lA)hi| c i(t)| 2 5 so that any deviation of 7 e //(t) from 70 is a signature of non- 
exponential decay. In addition, the eventual intersection J e tt(t) = 7o rules the transition 
from Zeno to anti-Zeno effects for repetitive measurements [9(. In order to determine the 
exact law for the survival probability P(t), one has to calculate the eigenfunctions of (1) and 
construct a suitable superposition of them corresponding, at t = 0, to a particle localized 
in the well |1), i.e. to c„{0) = 8 n> i. The tight-binding Hamiltonian (1) has a continuous 
spectrum of eigenfunctions [21( which can be calculated by separation of variables and cor- 
respond to c n (t) = u n {Q) exp[iQ(Q)t], where Q(Q) = 2cosQ is the dispersion curve of the 
tight-binding lattice band, — n < Q < n varies in the first Brillouin zone, and: 

u x = A(l + r)/(2cosQ), (5) 
u n = exp[— iQ{n — 2)] + rexp[iQ(n — 2)] {n > 2). 

In Eq.(5), r = r(Q) is the reflection coefficient for Bloch waves at the boundary of the 
semi-infinite lattice and reads explicitly: 

= A 2 - 2cosQexp(iQ) 
r lW A 2 -2cosQexp(-zQ)' 1 ' 

To study the decay process, we construct a superposition of the eigenstates, c n {t) = 
jZ n dQ F{Q)u n (Q) exj)[iQ(Q)t], where the spectrum F(Q) is determined by the initial con- 
ditions c n (0) = 5 n> \. Using an iterative procedure that will be described in detail elsewhere, 
one can show that the searched spectrum is given by F{Q) = — (27rA) _1 [A 2 exp(iQ) — 
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2cosQ]/[A 2 — 1 — exp(2i<5)]. Therefore the exact decay law for the occupation amplitude 
of site 1 1) is given by: 



Cl (t) 



1 

27 



f dQ exp(2^cosg)-l-4- Xp( 2 ' g) 
j-7r 1 + a 2 < 



(7) 



; exp(— 2iQ) ' 

where a = (1 — A 2 ) 1 / 2 . The short-time decay of |ci(t)| 2 is obviously parabolic; the long-time 
behavior of ci(t) can be calculated by use of the method of the stationary phase, yielding 
the oscillatory power-law decay 

1 1 - a 2 1 



Cl (t) 



cos(2t — 37r/4) as £ — > oo. 



(8) 



/tt(1 + a 2 ) 2 t 3 / 2 

In order to extract the exponential decay part from c±(t), after setting z = exp(iQ) it is 
worth rewriting Eq.(7) as an integral in the complex plane: 



cAt) = / dz exp 

w 2m L F 



it[z + - 

z 



(9) 



z(z 2 + a 2 Y 

where the contour a is the unit circle \z\ = 1. The integral (9) can be evaluated by use of 
the residue theorem. Note that, for A = 1 there is only one singularity at z — 0, and from 
residue theorem one obtains: 

Cl (t) = (l/t)Ji(2t), (10) 

which shows that, in the strong coupling regime, the decay greatly deviates from an ex- 
ponential law at any time scale. For A < 1, there are three singularities, at z = and 
z = ±ia, internal to the contour a. The residue associated to the singularity z = —ia yields 
an exponentially-decaying term, whereas the sum of residues at z = and z = ia yields a 
bounded function s(t), which can be written as a Neumann series. Precisely, one can write: 

ci(f) = VZ exp(- 7o t/2) + s(t), (11) 

where 7 is the natural decay rate as given by Gamow's theory [Eq.(4)], \/~Z = (a 2 + l) / (2a 2 ), 
and 

(20 



s(t) = J (2t) + ( 1 + \ 



1 ~ M2t) 
2,^ a 1 



E 



a 
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(12) 



is the correction to the exponential decay term. The decomposition (11) is meaningful in 
the weak coupling regime (A — > 0) since, in this limit, one can show that the contribution 
s(t) is small and of order ~ A 2 . The appearance of non-exponential features in the decay 
dynamics when approaching the strong coupling limit is clearly shown in Fig. 2, where the 



0.35 



200 



(a) 



(b) 



(c) 



10 15 20 25 
Time 



400 




200 



40 80 
Site Number n 



FIG. 2: Left: Behavior of the effective decay rate j e ft and amplitude |ci(t)| (insets) versus time. 
Right: Grey-scale image of |c n (t)|. In (a), A = 0.3; in (b), A = 0.5; in (c), A = 0.9. The horizontal 
dashed curves are the natural decay rate 70. 

numerically-computed behavior of the effective decay rate 7 e //(£) is shown for a few values 
of A, together with the temporal evolution of amplitudes |c n (t)|. The appearance of strong 
oscillations in the J e ff(t) curve when the coupling strength increases is a clear signature 
of an oscillatory decay dynamics which sets in even at intermediate time scales. Consider 
now the case of projective measurements of state |1) at time intervals t — r. In the weak 
coupling limit, where the decay deviates from an exponential law solely at short and long 
times, deceleration of the decay (Zeno effect) occurs for r < r*, where r* is the smallest 
root of the equation 7 e //(T*) = 70 for instance, for parameter values of Fig. 2(a) one 
has r* ~ 85 . In the strong coupling regime [Fig.2(c)], the decay is highly oscillatory, and 
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FIG. 3: Tunneling decay dynamics in a L = 50-mm-long semi-infinite waveguide array (left 
column) and corresponding discrete diffraction patterns (right column), (a) Weak coupling regime 
[a = 12 (im and oo = 16 /im, corresponding to A ~ 0.28]; (b) strong coupling regime [a = 12 [im 
and ao = 12.5 //m, corresponding to A ~ 0.86]. 



acceleration of the decay (anti-Zeno effect) may be observed for a value of r close to e.g. 
the first peak of 7 e //, where 7 e //('7") is larger than 70; for instance, for parameter values of 
Fig. 2(c) anti-Zeno effect may be observed for r ~ 2.34. In this case, repetitive observations 
correspond to suppression of the oscillatory tails in the decay process. 

Physical realizations of the tight-binding model (1) are provided by electron transport in 
a chain of tunneling-coupled semiconductor quantum wells [ijj] or by discrete diffraction of 
photons in a semi-infinite array of tunneling-coupled optical waveguides, where the tempo- 
ral variable t of the quantum problem is mapped into the spatial propagation coordinate z 
along the array [Fig. 1(b)]. Here we consider in detail the latter optical system since it shows 
several advantages: (i) Visualization of the tunneling dynamics is experimentally accessible 



18 



221 ]. and a quantitative measure of light decay can be done by e.g. NSOM techniques 



231 ] ; (ii) Preparation of the system on state |1) is simply realized by initial excitation of the 



boundary waveguide by a focused laser beam; (iii) Light diffraction experiments in waveg- 
uide arrays have successfully confirmed the reliability of the tight-binding model fl8| ; 
(iv) Transport of photons instead of charged particles (e.g. electrons) avoids the occurrence 
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FIG. 4: (a) Schematic of a waveguide array for the observation of the optical Zeno effect, (b) 
Numerically-computed behavior of mode amplitude |ci| trapped in waveguide |1) (solid curve) 
versus propagation distance in a L = 20-mm-long array for r = 4 mm, ao = 16 /im, and a = 12 (im. 
The dashed curve is the behavior corresponding to Fig. 3(a). (c) Grey-scale discrete diffraction 
pattern along the array. 

of dephasing or many-body effects, making waveguide-based optical structures an ideal labo- 
ratory for the observation of several analogs of coherent quantum dynamical effects (see, e.g. 



22], Landau-Zener tunneling 



). Beautiful optical analogs of Bloch oscillations |14. Lsl . 
22j, adiabatic stabilization of atoms in strong fields [24]], and coherent control of quantum 
tunneling 25j, have been indeed reported in recent optical experiments. 
Light propagation in the waveguide array is described by Eq.(2) in which the temporal vari- 
able t is replaced by the spatial propagation coordinate z, % = \/(2ir) is the reduced wave- 
length of photons, m = n s is the refractive index of the array substrate, V(x) ~ n s — n(x), 



and n(x) is the array refractive index profile (see, e.g., [24j, [25|]). As an example, Fig. 3 
shows the discrete diffraction patterns and corresponding behavior of light trapped in waveg- 
uide |1) as obtained by a numerical analysis of Eq.(2) using a standard beam propagation 
method with absorbing boundary conditions 26j| ; initial condition corresponds to excitation 
of waveguide |1) in its fundamental mode, i.e. ip(x,0) = f(x). The refractive index profile 
of the semi-infinite array used in the simulations is plotted in Fig. 1(c) for parameter values 
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which typically apply to lithium-niobate waveguides |24j. Note that, as A is increased, non- 
exponential features are clearly visible. However, as compared to the tight-binding results, 
the peaked structure of J e ff(t) obtained from the continuous model (2) is smoothed [compare 
e.g. Fig.2(c) and Fig.3(b)]. In order to reproduce the optical analog of the quantum Zeno 
effect in the waveguide system, one can adopt the array configuration shown in Fig.4(a), in 
which a straight waveguide |1) is periodically coupled, at equally-spaced distances z = r, to 
semi-infinite arrays of finite length r placed on alternating sides of the waveguide. At each 
section where the lateral arrays end, light trapped in the interrupted waveguides is scattered 
out and solely a negligible fraction of it will be re-coupled into the waveguides at the next 
section of the array. Therefore, at planes z = t, 2t, 3t, ... one can assume, at first approx- 
imation, that a collapse of the state ip(x, z) into the fundamental mode <p(x) of waveguide 
|1) occurs, thus simulating the 'wavepacket collapse' of an ideal quantum measurement. An 
example of deceleration of the decay via tunneling in the alternating array, analogous to the 
quantum Zeno effect, is shown in Figs. 4(b) and (c). 

In conclusion, an exactly-solvable model for the tunneling escape dynamics of a metastable 
state coupled to a tight-binding lattice has been presented, and its optical realization - 
including an optical analog of the quantum Zeno effect- has been proposed in an array of 
tunneling-coupled optical waveguides. 
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